Qubit purification speed-up for three complementary continuous measurements 
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We consider qubit purification under simultaneous continuous measurement of the three 

non-commuting qubit operators a x , & v , a z . The purification dynamics is quantified by (i) the 

Oh' average purification rate, and (ii) the mean time of reaching given level of purity, 1— e. Under 

ideal measurements (detector efficiency 77 = 1), we show in the first case an asymptotic mean 

purification speed-up of 4 as compared to a standard (classical) single-detector measurement. 

However by the second measure — the mean time of first passage T(e) of the purity — the 

corresponding speed-up is only 2. We explain these speed-ups using the isotropy of the qubit 

evolution that provides an equivalence between the original measurement directions and three 

0^ ■ simultaneous measurements, one with an axis aligned along the Bloch vector and the other 

00 

£C) ■ with axes in the two complementary directions. For inefficient detectors, r/ = 1 — 5 < 1 

the mean time of first passage T{5, e) increases since qubit purification competes with an 

o ; 

£NJ , isotropic qubit dephasing. In the asymptotic high-purity limit (e, 5 <C 1) we show that 

1 _ 

the increase possesses a scaling behavior. AT(5,e) is a function only of the ratio 5/e. The 

increase AT (5/e) is linear for small argument but becomes exponential ~ exp((5/2e) for 5/e 

X ■ 
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PACS numbers: quantum measurement, quantum control, quantum feedback, qubit, purification 



I. INTRODUCTION 

Pure states are an important resource in quantum computation and communication algorithms 
[lljQ]- While state purification is possible via cooling, this may be impractical for several reasons, 
including long relaxation times, degenerate ground states, and the presence of several dephas- 
ing mechanisms. A different purification process is possible via continuous measurement when 
the (available to the observer) quantum state will purify continuously according to the detector 
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measurement result. Here, the speed at which one can purify the state is set by the detector mea- 
surement rate, and the final purity will depend on how close the detector is to a quantum-limited 
(100% efficient) one. This may become an important tool since continuous measurements are also 
at the heart of various quantum control applications [2J1, including quantum state stabilization via 

nn n 

quantum feedback [3H5J, preparation of entangled states |aj, and for continuous error corrections 

A 

In recent years several groups have suggested rapid purification protocols based on continuous 
measurement and Hamiltonian feedback [8l4l2|. which makes it possible to considerably speed-up 
purification. For a single qubit, the problem was first analyzed by Jacobs 8], who recognized that 
the rate of average purification can be enhanced using unitary transformation at each measurement 
time step, so as to make the state always orthogonal (in the Bloch sphere sense) to the detector's 
measurement basis. This feedback algorithm (which has been rigorously shown to be optimal 
and rederived in Ref. [l4]) produces a factor of 2 speed-up in the high-purity limit, when the 
state approaches the Bloch sphere surface. The speed-up is by comparison with a no-feedback 
measurement, which can be completely understood classically (see below). This speed-up is a 
quantum mechanical effect since it is only possible if the system can exist in superposition of the 
eigenstates of the measured observable. 

Jacobs' protocol maximizes the average purification rate in the sense that it minimizes the 
time r it takes the average purity {p) to reach a certain level 1 — e. A different, but equally well 
motivated, aim was suggested by Wiseman and Ralph [llj, that of minimizing the mean time (T) 



for the purity p to attain a certain level. It was argued (and later rigorously confirmed [13l . Il4l ]) 
that the optimal protocol in this case is to rotate the state so that at each measurement time step 
it is aligned with the detector (i.e., making the state density matrix diagonal in the basis of the 
detector observable). This is the classical protocol referred to above, and in the absence of any 
other dynamics, requires no feedback to realize. 

In this paper we study a recently suggested purification protocol based on monitoring of the 
qubit state via simultaneous continuous measurement of the three non-commuting qubit operators 



151 ] a x , a y , a z (FigQ]). (A related protocol based on simultaneous measurement of these observables 
was first considered in Ref. In contrast to the above purification protocols there is no need 
to perform Hamiltonian feedback. In this sense it is an "open loop" protocol, similar to the 
"random unitary control" purification protocols introduced in Ref. [13] for spin-j systems. Though 



operationally different, the two protocols can be shown to be equivalent in the special (j = 1/2) 
case, which is another consequence of the isotropy of the qubit evolution discussed below. 
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FIG. 1: A qubit measured by three orthogonal detectors. 



Considering the rate of average purification under ideal measurements (with detector efficiency 
rj = 1) the purification speed-up is 4, (this speed-up is implicit in the results of Ref. (3]). Here 
the comparison is with a single measurement of the same strength with no feedback (the speed-up 



of 4/3 quoted in Ref. [1 71 ] comes from keeping the total measurement strength the same in the 
two protocols). However, we show here that by the alternate metric — the mean time (T) of first 
passage — the speed-up is only 2. Both of these purification speed-ups can be understood via 
the isotropy of the qubit evolution [jjj in the Bloch space that allows one to represent the three 
detectors by another equivalent detector triad, <5v, <5y, <5y, measuring at each time moment along 
the state and in the complementary directions. 

We also study the purification dynamics for inefficient detectors, rj = 1 — 5 < 1. In the 
asymptotic high-purity limit (e, 5 <C 1) the increase AT (5, e) of the mean time of first passage (T) 
shows a scaling behavior (i.e. it depends only on the ratio 5/e) as seen from numerical calculations. 
The increase AT (5/e) is linear for relatively small inefficiency but grows exponentially, reaching 
~ exp(5/2e) for relatively large inefficiency. 



II. FURTHER BACKGROUND: SINGLE DETECTOR PURIFICATION PROTOCOLS 



A. Qubit impurity from single-detector measurement 



For a quantum-limited detector, in the limit of infinite detector bandwidth, the evolution of 
the state p(t) of a quantum system due to weak continuous measurements of a variable X can be 
described by the stochastic master equation (SME) Q]: 



dp = ^ dtV[X] p + ^dW(t)H[X] p. 



(1) 



4 



Here V[A)p = ApA^ - \{A^ Ap + pA^A), H[A]p = Ap + P A* - Tr[(A^ + A)p]p, and dW is the 
increment of a Wiener noise process. The internal Hamiltonian evolution of the system is either 
neglected (in case of a large measurement strength To), or eliminated (see, e.g. Ref. Q). The 
measurement rate Tq sets the rate at which information (about the final state) is extracted, and 
thus the rate at which the system is projected onto a single eigenstate of X fla . [3] (provided the 
spectrum of X is non-degenerate). We will specialize on the measurement of a qubit and choose 
first X = a z , although some results below are valid for a general X. 

The measurement record in a small time interval [t, t + dt) is given by 



dR(t) = (X) dt + dW(t)/y/2To, 



(2) 



where {X) = Tr[X p(t)], and dW(t) is the same realization of the Wiener increment that appears 
in Eq. (fTJ, thus providing explicit dependence of the state evolution on the specific measurement 
result. The evolution gene rated via Eq. ([I]) corresponds to infinitesimal measurement operators 
(see, e.g. Refs. [2, 
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23l|). Since the same observable X is measured at each time step, these 



operators commute, allowing one to integrate Eq. (fTJ exactly. Thus, for a finite time interval [0, r) 
the final state will depend only on the time-averaged measurement record: 



p( T ) = - I dt'dR{t'). 
T Jo 



(3) 



Introducing the measurement operator 



M, 



Tor 



TT 



1/4 



exp [—((J, — ct z 



,2 T 0T, 



one obtains, using the method of un-normalized density matrices (see, e.g. Refs. 

M^ z p{0)Ml z 



20 23]): 



(4) 



Here the probability distribution for p{j) is given by P{p) = TrfpM^Mu J. 

In the basis of the observable X, the above result can be written as a quantum 

nn 

Bayesian filter [181 . 1241 ]: the update of the diagonal density matrix elements pa will look ex- 



actly as a classical Bayesian update of a "probability distri 
with Gaussian likelihoods P(p\i); this was termed in Ref. 



bution": Pu (t) = Pii (0)P(p\i)/P(p), 
181 ] the quantum- classical correspon- 



dence principle. The outcomes p will be Gaussian distributed around the eigenvalues of X, {xi}: 



P(p\i) = ^/Tot/tt exp[— (/i — Xi) 2 ToT] with variance Var = l/2For, and expected outcome probabil- 
ity P(p) = PuP{lA^)- The update of the non-diagonal matrix elements in ^ will be according 



to the rule: pij{r) = pij (0) yj pa (r)pjj (r) / pa (0) pjj (0) provided the measurement is preformed by a 



quantum-limited detector; this implies that a pure state will remain pure From the above it is 
clear that the measurement time T me as to distinguish approxima tely two eigenvalues, set 
by Var| Tmcas = (xi - xjf/A. For a qubit this gives T meas = l/2r [18|. 



For further use we also write down the quantum state evolution in terms of the Bloch vector 
components of the state, x = 2Hepi2, y = 2Im/9i2, z = p\\ — P22- From Eq. ([I]) for the measurement 
of a z it follows: 

dz = (1 - z 2 ) a/2T7 dW 



o 
p 

dx = x — zx \/ZTq dW 

V 

dy = -^y-zy^/2T^dW. (5) 

n 

In these equations we also included the effect of detector non- ideality (inefficiency) r], leading to a 
pure dephasing : In the ensemble averaged equations (since Eqs.([5]) are in the ltd form, averaging 
means just to drop the noise term, and corresponds to ignoring the detector results), Tq is the 
decoherence rate due to an ideal detector of measurement rate Tq, and rj = Fo/(ro + 7) < 1 is the 
detector efficiency (ideality), defined as the ratio of Tq to the total decoherence r = Tq + 7- For 
a single detector the simplest model to describe the extra decoherence 7 is to consider a second 
independent detector "in parallel" (i.e., measuring the same a z variable) by adding terms similar 
to Eq. (P) with Tq replaced by 7, and subsequently averaging over that detector output. Thus, 
the density matrix available for an observer who takes into account only the results of the first 
detector will be described by Eqs.(j5]). 

B. Single-detector purification protocols 

We consider first purification protocols via single-detector measurement with an ideal (quantum- 
limited) detector. 

1. Purification without feedback 

For a single detector measurement without feedback, the state evolution in the detector basis is 
essentially classical, Eq. That is, if the density matrix begins diagonal in the measurement basis 
(as it will be if it is a completely mixed state), it remains so and there is no way to distinguish 



{pn(t), p22{t)} from a classical probability distribution [18|, |2a|. For continuous measurement 



one can consider the purity or entropy of the monitored state, and in this particular case the 
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von Neumann entropy of the state coincides with the Shannon entropy (see, e.g. Refs. 
SVn = — Tr[p In p] = — ipxi l n Pli + P22 hi ^22)- In what follows we will consider the so called linear 
entropy (see, e.g. 26]). s = 1 — Trp 2 = 1 — p which is a monotonic function of S v n- Here 



'vN- 

p = Tr/r = 2^ is the purity expressed through the Bloch components of the state. 

The corresponding equation for the purity (at i] = 1) follows from Eq. ([5]): 



dp = 2T [(1 - p)(l - 2 2 )] + 2 V / 2T^z(l - p) dW. 



(6) 



Since purity p, like entropy, is invariant under unitary transformations, without loss of generality 
one can say that measurement "parallel to" (in the same basis as) the state corresponds to x = y = 
and 



dp = T [1 - z 2 ] 2 dt + y/2T Q z[l - z z ] dW. 



(7) 



so that the average change of purity is (dp}\\ = Tq (1 — z 2 ) 2 dt. We note that the same result holds 
for non-ideal measurement, rj < 1, since the detector's non-ideality affects only the evolution of the 
non-diagonal density matrix elements; this means that a non-ideal detector will purify the state 
if the measurement is along the state. It is clear that Tq plays the role of a maximal classical 
purification (information acquisition) rate that happen when z = p\\ — P22 = 0, i.e. when the two 



outcomes are equally likely (see, e.g., Ref. 



281]). By approaching z — > ±1, (dp) 11 — > since "little 



information" remains to be extracted to clarify that collapse has happened 



2. Jacobs feedback purification protocol 

From the discussion above, it is intuitively clear how the Jacobs' enhanced purification protocol 
works. Given the state p(t), one should continuously adapt the measurement basis [or, equivalently, 
rotate the state to some p'(t)\ so that the detector would perform measurement in a complementary 
direction with respect to the eigenbasis of the rotated state (i.e., measuring perpendicular to the 
state, in the Bloch picture). In the detector basis the rotated state again possesses equally likely 



281 ]. this procedure maximizes the average 



outcomes, with p'n(t) = p^ify- As expected classically 
purification rate since p = p' under rotation, while z — > z' = in Eq. (|6|). However, the possibility 
to perform coherent rotations of the density matrix is of course a quantum mechanical effect. The 
procedure also makes the purity evolution deterministic jsj], i.e. the noise term in Eq. ([6]) is zeroed 
so that 



dp = 2r (1 - p)dt 



(8) 
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and {p) = p under this protocol. From this one can evaluate the time tj_ when the average purity 
[p) reaches a given level 1 — e: 

It should be noted that an attempt to evaluate the analogous time tm in the case of the classical 
measurement along the state, using {dp)u from Eq. ([7]), will lead to a wrong scaling of ~ 1/e. 
The reason is that (dp) ^ dip) in general, and the evolving purity distribution V(p,t) is generally 
different from <5-function [10(]. The correct evaluation of tii is to calculate the average purity {p)(t) 
by taking into account the exact solution, Eq. (jl]), of the stochastic evolution equations ([!]). In the 



high-purity limit this leads to 



-Lure" 1 (10) 
1 o 



This is exactly twice as long as the time in Eq. Q, which establishes the speed-up of 2 for the 
Jacobs protocol. We note however, that, unlike the case of parallel measurement, the perpendicular 
measurement will not purify to a completely pure state if the measurement is inefficient, as will 
be explored in Sec. IIIICI The reason is that the excess back-action in an inefficient perpendicular 
measurements results in decay in the coherences of the state. 

3. Wiseman and Ralph purification protocol 

Instead of considering the time r at which the average purity (p) reaches a certain level 1 — e, 
one can also consider the average time (T) for a system to attain that purity level, p(T) = 1 — e 

inn □ 

Id . |12l . Il3l ] . It was noted in Ref . [10J that in many circumstances the time (T) is a more useful 
quantity. The reason is that the time T at which p(T) = 1 — e, has a well-behaved statistics [ic3 ]. 
in contrast with p(t) which has extremely long tails at relatively small values of p. That is, the 
averaged (p) is strongly influenced by the rare cases that are slow to purify. Because of this there 
is a substantial disagreement between r and (T) for a qubit. It was shown in Ref. [hJ, however, 
that good agreement is found between (T) and T log , defined as the time required for (ln[l — p(t)]) 
to reach the certain level hie. This is because taking the logarithm de-emphasizes the tails, and 



indeed for a qubit ln(l — p) has near-normal distribution [1C|. 

Therefore, we consider the stochastic equation for the logarithm of the linear entropy s = 1—p. 
It follows from Eq. ([6]) that: 

dins = -2T {2s + x 2 + 2z 2 }dt + 2y/2T^zdW (11) 
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(Here, without loss of generality we have put y = 0, since a single admeasurement keeps the 



state of the qubit in a fixed meridional plane). In the Wiseman- Ralph feedback protocol Kj one 
keeps the monitored state along the detector z-axis (so x 2 = 0, while z 2 is maximized at 1 — 2s). 
Thus, the single detector purification is maximized and in the high-purity limit, s ~ 0, we obtain 
— ~ — 4ro- Using this, the average time (of first passage) (T) is evaluated as (T) ~ T log and 
in the high-purity limit, E<1 

T^ = -Ll n£ -\ (12) 
II 41 o 

which is half the size of t±, and 4 times shorter than ry . In what follows we will use these results to 
understand the purification via three simultaneous complementary measurements without feedback. 

III. CONTINUOUS MEASUREMENT OF THREE COMPLEMENTARY QUBIT 

VARIABLES 

We consider continuous measurement of the qubit complementary observables a x , & y , and a z by 
three independent linear detectors (see Fig. [I]). In principle, given a <7 2 -detector, the measurement 
of, say, & x can be implemented operationally by performing fast unitary rotation, U x (t), of the 
state towards the x-axis "at the beginning" of the measurement interval dt, then a continuous 
measurement with the z-detector and finally, a backwards transformation, JJ^jt + dt) "at the 



171 ] . We also mention 



end" , in any infinitesimal measurement time step dt, as discussed in Ref. 
that a simultaneous measurement of the three observables is possible to implement, at least in 
principle, in a quantum optics setup j^] (e.g., using a single atom in an optical cavity field). 

In the case of detectors measuring the qubit in the bases of non-commuting observables, it is 
impossible to use the quantum Bayes rule ^ for finite times because the measurement back-actions, 
M^x, May, Af„ z , do not commute with each other. The measurement back-actions commute 
only for small measurement time intervals dt so that one should apply the POVM update in its 
differential form Eq. ([1]) and then sum up the contributions to the qubit evolution. For the case of 
mutually unbiased observables & x , & y , and a z it is convenient to use the Bloch vector components, 
r = (x,y,z) = {Tr\pak],k = x,y,z} determined by the time-dependent qubit density matrix p. 
From Eq. ([5]) for the influence due to (^-measurement, the influences due to a x (cr ?/ )-measurement 
can be obtained by cyclic replacements z — > x — > y, bearing in mind that the variables x(y) should 
play the same role as z in a <7 2 -measurement. In this way we obtain the following evolution equations 
in Ito form (see Ref. [15I ]) that take into account the three measurement records dRk(t), k = x,y,z 
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as in Eq. ([2]), and correspondingly introduce three independent Wiener processes, dW% dW^ = 5ik dt: 

dx = -(T y + T z ) x + (1 - x 2 ) V / 2TV: dW x 

-xy^/W^ y dW y -xz^/W^ z dW z (13) 
dy = -{T z + T x ) y + (1 - y 2 ) ^/W^ y dW y 

-yz^2l\ z dW.-xy^W^ dW x (14) 
dz = -(T x +T y )z + (1- z 2 ) ^/W^ dW z 

-x z v / 2To~^ dW x -yz ^/W^ y dW y (15) 

Here Tf, = Jk+ ^o,k are the total decoherence rates for each detector, including individual dephas- 
ings 7^ and measurement rates To^,k = x,y,z. Note that a pure dephasing, say in the x-basis, 
causes the decay of the z- and y-components. Eqs. (j!4[) and (I15p are simply obtained by cyclic 
permutation of the variables in Eq. (fl3j) . 



A. Qubit evolution with identical detectors 

In what follows we consider the simplest (but still rich) case of three identical detectors: Tq^ 
To and 7fc = 7 > 0. Then the qubit evolution (|13[) — ()15j) can be rewritten in a vector form as 



dr = -2T rdt + {dW(t) (1 - r 2 ) - [r x [r x dW(t)]\ } , (16) 

where dW(t) = {dW x ,dW y ,dW z } is the vector of Wiener increments corresponding to the vector 
of results dR = {dR x ,dR y ,dR z }, similar to Eq. ([2]). The evolutio n (1161) is invariant under arbitrary 
rotations (see Fig. [2]) of the coordinate system in Bloch space 3CJ] (as is the ensemble averaged 
evolution: r = — 2Tr). Hence while measurement of only a single observable "attracts" the 
q ubit state to one of the c_ d ,„ g e, g en TO t„ rs , the sh-taneous _e„t of »., » , 
and a z leads to no preferable direction in the Bloch space. In Ref. [15J this isotropy of the qubit 
evolution was shown to lead to a locally isotropic Brownian diffusion of the direction of the Bloch 
vector. In particular, for pure states under ideal measurements, the point on the Bloch sphere 
diffuses isotropically 



311 ] with a diffusion coefficient Tq [15j . The evolution isotropy was used then 
to construct simple quantum state estimations [3]. In what follows, we will use the qubit evolution 
isotropy to understand the qubit purification dynamics. 



dR « dW 




FIG. 2: Bloch equations invariance under rotations for three complementary measurements with identical 
detectors. The evolution is determined by the relative orientation of the Bloch vector and the vector of 
results 151. 



B. Purification dynamics with three complementary measurements 

Eq. (|16p shows that the measurement term component along the radius is \/2Tq (dW-e r ) (1 — r 2 ) 
and vanishes on the sphere. The decreasing r-dependent coefficient suggests that r = 1 is an 
attractor of the random evolution. Thus, the qubit state will purify in the ideal case (7 = 0) for 
any realization of the measurement process. 

To establish further the purification dynamics, we start with Eq. (|16p and transform it to polar 
coordinates, r, 6, 4>. The equation for the radius decouples from the other two and reads [15I ]: 



dr = 2T (1/r - r/r])dt + V2To(l - r 2 ) dW r 



(17) 



where dW r = dW ■ e r is the projection of the noise vector onto the Bloch vector direction, 
e r = (sin 9 cos (p, sin 9 simp, cos 9), and the normalization of the noise is dW r dW r = dt. 

As seen from Eq. (|17p . if r < r] 1 ^ 2 , then on average r > 0, which means state purification for ideal 
measurements, 77 = 1. Since the equation for r is singular at the origin, it will be more convenient 



to consider further the dynamics of the purity p = Tip 2 
a totally mixed state corresponds to p = 1/2. We find 



(l+r 2 ) 



In this choice 1/2 < p < 1, i.e., 



dp = 2T {1 - (2p- 1)/t/ + 2(1 -p) 2 }dt + 2y/2T~ (l -p) v / 2p- ldW r . 



(18) 



For an ideal measurement, i] = 1, and starting from a non-pure initial state, the qubit will purify 
(p — > 1) on a time scale of the order of r mea s ~ ^0 ■ For a non-ideal measurement, 77 < 1, purity 
will continue to fluctuate around a stationary average value, {p) s t < 1- 
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1. Fokker- Planck equation for purity 
To quantify the purity dynamics and asymptotic purity distribution we first consider the Fokker- 



Planck equation (FPE) corresponding to Eq. (118|) . Using the standard coefficients (see Ref. 321]): 



A(p) = 2T [1 - (2p - l)/r? + 2(1 - p) 2 ] (19) 
B(jp) = ZT Q {2p-l){l-p) 2 (20) 



the FPE reads 



and the initial distribution is taken to be V(p,0) = 5(p — po). At t T meas the purity reaches a 
stationary distribution [lj| 

(2p-l)(l-77)- 



.! V2^T 
(l-p)S 



^st(P, ? ?) = ^ Ti 1\3 eX P 



(21) 



2(1 -p)»7 

where iV is the normalization. For rj — >■ 1, the stationary distribution V s t{p,r]) approaches the 
5- function at p = 1. 

The purification dynamics can be approximated using the ensemble-averaged purification rate 
[s], [l(3] obtained from the above ltd equation (fT8|) for purity: 



(dp) = 2T (1 - (2p - l)/rj + 2(1 - pf) dt. (22) 

Now consider a naive approach to integrating Eq. (|22p . in which we replace p by (p) everywhere 
it appears on the right-hand-side: 

(dp) = 2T [1 - (2(p) - I) /n + 2(1 - (p)) 2 } dt. (23) 

Contrary to the single-detector measurement, it can be shown, by numerically solving the FPE, 
that the evolution of the average purity (p) obtained from the naive Eq. (|23p is very close to the 
true average (p(t)) FP = J 1/2 pV{p,t)dp, obtained from FPE. (The reasons for this difference will 
be discussed below.) Numerically, the approximation is best in the high purity limit ((p) = 1 — (s), 
(s) 1) and for almost ideal detectors, rj ~ 1. In the high-purity limit we can discard terms 
of order (s) 2 (i.e. the final term) so that Eq. (|23|) becomes a simple linear equation. For ideal 
measurements the time when the average purity reaches the level of 1 — e, e <C 1 is thus 
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which is four times shorter than the standard time rn of classical measurements, Eq. (jlOp . This 
means with three orthogonal measurements we obtain 4 times speed-up as compared to 2 times 
speed-up of the Jacobs (single measurement) purification protocol. 

This result can be understood as follows. By the isotropy of the qubit evolution under three 
complementary measurements, in each time moment t one can represent the measurements in the 
x, y, z directions by measurements in an equivalent triple of directions x', y' , z' (in the Bloch 
space, Fig. [2]), chosen so that z' is parallel to the state (i.e., in the basis of <5v the density matrix 
p(t) is diagonal), while directions x', y' are perpendicular to the state; of course these chosen 
directions must change in time according to the evolution of the state, p(t). The measurements in 
the x' , y' directions in the time interval [t, t + dt) are termed "good" measurements as they are 
precisely the unbiassed measurements of Jacobs that maximize the average purification rate, and 
the ^'-measurement along the state would be termed "bad" by obvious reason. This observation 
could be confirmed by rewriting Eq. (|22p for r\ = 1 in the following way: 

(dp) = 2T ((1 - p) + (1 - p) + 2(1 - pf)dt, (25) 

where the first two terms correspond to the deterministic change, dp, Eq. ([8]), and the last term 
coincides with the average change of (dp) ii due to measurement along the state, Eq. (|7|). In the high 
purity limit, the last term is suppressed and the two "good" measurements contribute each a speed- 
up of 2 that amount to a total speed-up 35[ of 4. Moreover, this makes clear why the naive approach 
(replacing (p 2 ) by (p) 2 ) to calculating the mean purity works for the isotropic measurement: Jacob's 
protocol gives deterministic growth of the purity, so that (p 2 ) = (p) 2 . For isotropic measurement 
the purification in the mean is dominated by the perpendicular measurements, as just shown. 
Therefore the purification is approximately deterministic, and becomes more so the purer the state 
becomes. By contrast, with a single measurement and no control (the parallel measurement case), 
the purification is far from deterministic. In particular the long tail of low purities means it is 
impossible to obtain accurate results by replacing (p 2 ) by (p) 2 1£\. 

The "splitting" of the measurements into "good" and "bad" is also applicable when the goal 
of minimizing the mean time (T) is considered. As in the single-detector case JjJ, we take the 
log-entropy evolution, Ins, s = 1 — p; similarly to that case Ins should have more symmetric 
distribution and therefore its average rate of change would correspond to the mean time of first 
passage (MTFP). Using ltd equation for three ideal detections, the average change of log-entropy 
in the high-purity limit (s — > 0) reads: 

(dins) ~ -8T dt = -(4 + 2 + 2)r (ft, (26) 
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i.e., the three detector measurements now "splits" into one "good" measurement (V), directed 
along the state, which gives the rate 4ro as in Eq. (I12p [Wiseman and Ralph protocol], and two 
"bad" measurements (V, y'), directed perpendicular to the state, that give one-half of this rate 
each, see Eq. ([9]). Correspondingly, the mean time (T) ~ T log for the average log-entropy (Ins) to 
reach hie is given by 



This is two times shorter than the mean time for a single-detector no-control measurement, Tj, og , 
Eq. (|12p . Therefore, in comparison to the parallel single measurement case, three complementary 
measurements give a speed-up of 2 in terms of the mean time to attain a given purity. We will 
verify this result by explicitly calculating the MTFP, in the following subsection where we consider 
non-ideal detectors. 



It is now interesting to ask the question: Given non-ideal detectors with rj = 1 — 5, what level of 
purity, p = 1 — e can be reached and what time is needed? We now consider in detail this question, 
with emphasis on the high-purity and high-efficiency limit, when e <C 1 and 5 <C 1. The answer to 
the above question will depend on the goal examined under purification. 



Since the average purity {p(t))pp described by the FPE is numerically close to that from the 
naive ensemble-average evolution Eq. (|23p . (p(t)) (see also Ref.flj]]), we will use the latter in our 
analysis. In particular, the stationary value (p) s t is close to the true stationary value, obtained 
from the stationary distribution, Eq. (|2ip . and reads: 




(27) 



C. Purification dynamics for non-ideal detectors 



1. The goal of having the average purity reach the level (1 — e) 




(28) 



In the high-ideality limit it gives {p) s t — 1 — f • Therefore, (p) can reach a purity level (1 — e) only 



if 




(29) 



14 



The time for (p) to reach 1 — e can be calculated from Eq. (|23p . and in the high-purity limit it 
gives: 

T iso (5,e) ~ — 

41 o 

Here we note that the purification rate can be again understood by splitting into "good" and "bad" 
measurements. Indeed, the single detector average purification for rj < 1 can be calculated from 
Eq. ©: 

(dp) = 2T /(l - p)(l - z 2 ) + Ul - ±) (2p - 1 - z 2 )^ dt, (31) 

and one can observe that (dp) is maximized for a measurement in perpendicular direction (z 2 = 0) 
for not too small rj (e.g. even for a totally mixed state, p = 1/2, this happens for rj > 1/2). 
Thus, the first term in Eq. (|22p comes exactly from two "good" measurements (perpendicular to 
the state), 2(dp)± = 1 — (2p — l)/r/, while the second term comes from the "bad" measurement 
(parallel to the state) as in Eq. ([7|), and is suppressed in the high-purity limit. Therefore, the 
asymptotic result for Ti so (S, e) comes entirely from the two good measurements, as in the ideal 
case. The time increase in Eq. (130p with respect to the ideal measurement, Eq. (I24j) is linear in 
the ratio 5/e for 5 <C e and diverges logarithmically to oo as the impurity level e approaches the 
bound, Eq. ([29]) . To reach impurity level e < | is simply impossible. 



hrl-lnfl-A 

2s \ 2e 



(30) 



2. The goal of having a certain mean time (T) at which the purity p(T) reaches the level (1 — e) 



We now proceed with the investigation of the mean time (T) for t 
for three complementary non- ideal measurements. Following Ref. 
for the log-entropy: 



le purity to reach p(T) = 1 — e 



3 



we write the Ito equation 



dins | iso = -2r 1 2 - 2s + ~ + ^1 - i j dt + 2 y^Vl - 2s dW r . (32) 

Similar to the ideal case, the average change of log-entropy can be represented as a sum of the 
measurement along the state and two measurements perpendicular to the state. Indeed, the change 
dins for a single non- ideal detector can be written as: 

dlns = -2r ^2s + x 2 + 2z 2 + (l - |-| dt + 2 v / 2T^z dW. (33) 

We note that for given entropy s = 1 — p, the average change (din s) is maximized again for x 2 = 
and z 2 = 1— 2s, and it coincides with the ideal case: (dlns)n = — 2T$dt {2— 2s}, see Eq. (1121) . On the 
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other hand the change due to measurement in the complementary directions (z = 0) is deterministic 
and read: (d\ns}± = <ilnsj_ = — 2Todt [- + (1 — -)/2s], so that (dins) |i so = (^ ms )|| + 2 (cHns)_i_. 
By considering highly ideal detectors, rj ~ 1, so that the last, singular, term in dlns± can be 
neglected, one can reproduce the mean time (7 1 ) by integration of the equation for (dins) |iso* 
However, in the high-purity limit this is not possible: naive integration of log-entropy, with the 
singular term, oc (1 — r)/s, included, will lead to a wrong result. Since in the high-purity limit 
6 <C 1, e <C 1, this term is of the order of |, it is clear that the time of reaching a certain purity 
level will be essentially affected when | > 1. 

In what follows, we investigate the exact solution for the MTFP in case of inefficient detectors. 
We consider an initially completely mixed state, and denote the MTFP as T(8,e). The FPE 
constructed above, based on the stochastic ltd Eq. (|18p is used, where in the high-purity limit, 
we approximate the coefficient, Eq. (|20p . as: B(p) ~ £?hp(p) = 8Fq (1 — p) 2 , while keeping A(p) 
unchanged. The MTFP solution for one reflective boundary (at p = 1/2) and one absorptive (at 



p = 1 — e) is written as (see, e.g. Ref. 



Jpo ^hp(2/, 1-5) J 1/2 Bnp(z) 
where po is the initial purity as before; in what follows we consider only po = 1/2. The function 
that enters the solution is ^hp(^) = exp Jy 2 dx' J~7^7) and in our case one obtains 

exp (x — i) fi^s(x) , . 

[2(1 - x)]^) 

where we denoted fis(x) = 1 - (1 _ x) <5 (1 _ [5) ■ 

For ideal measurements, 5 = 0, one can integrate Eq. (|34|) to obtain analytically the result: 

f(0,e)^-L (Ins" 1 - 1.35), (36) 
81 o 



in agreement 



33J] with the straightforward integration of the ensemble- averaged log-entropy equa- 
tion, [see Eqs. (l26p and (|27p above]. For 5 > analytical integration is not as obvious: naive 
expansion of the integrand in series in 5 leads to a sum of contributions with increasing singulari- 
ties. We present instead numerical calculations that show the scaling behavior in the high-purity 
limit. 

The main observation is that in the high-purity limit the time of purification can be represented 

as 

r(*, e )~r(o,e) + Af^y (37) 
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FIG. 3: The mean time increase, AT(ae,e), for small a < 5, and £« 1; Scaling is established for e < 10 
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FIG. 4: The approximate linear growth of the logarithm of mean time increase, ln[AT(ae, e)] ~ C\a 1 for 
large a 3> 1, and e< 1; Scaling is for e < 10 -4 . 

That is, the time increase AT depends only on the ratio a = |. It can be shown that for a <C 1, 
AT (a) ~ (l/8r )§, while for a > 1, AT (a) » (l/8r )exp [Ci(a) a], with 0.25 < Ci(a) < 0.5, 
approaching 0.5 for large a. 

On Figs. [3] and H] we present evidence of scaling. Fig. 3 shows AT (a) = T(ae, e) — T(0, e) at 
fixed e and for a < 5. The curve for e = 10~ 3 is close to the curves for e = 10 -4 , 10~ 5 , 10~ 6 that 
practically coincide on the figure; the slope at a = is 1/6. Fig. 4 shows the scaling for large a. The 
ln[AT(a)], eventually approaches a single curve for e«1 For a > 20 it looks like a linear growth, 
i.e. ln[AT(a)] ~ C\a. More precisely, by numerical calculation of the derivative d\a\AT(a)\/da 
(not shown) one can see that it is not a constant but varies somewhat as discussed above. For 
small enough e (a > 2000), C\ approaches 1/2. Thus, the numerical calculations support Eq. (I37|) . 
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IV. CONCLUSION 

In this paper we investigate in detail a recently proposed quantum state purification protocol 
based on simultaneous continuous measurement of the three complementary observables for a 
qubit: a x , a y , a z . Contrary to analogous single-detector purification protocols, that are based on 



complementarity Hid, |34|, and which require feedback control, here there is no need of introducing 



quantum feedback. However, interestingly enough, the purification dynamics of our protocol can be 
understood via the above mentioned protocols, using the established isotropy of the qubit evolution 
under three complementary measurements by identical detectors [if]]. 

For ideal (i.e. quantum- limited, or efficient) measurements, our main results are the observation 
of different factors for the speed-up (relative to a single measurement in the eigenbasis of the qubit 
state), of 4 and 2, depending on how the purification speed is quantified. In the first case it is 
in terms of the time when the average system purity, (p), reaches certain purity level 1 — e; in 
the second case it is in terms of the mean time (T) at which the purity p(T) first attains the set 
level of 1 — e. Both of these speed-ups can be understood via the possibility to "split" the three 
detector measurement at each measurement time step into three equivalent measurements — one 
parallel to the state (in the Bloch sphere sense) and the remaining two perpendicular to the state 
- which is one of the consequences of isotropy of the qubit evolution in the Bloch space. For the 
first measure of purification time these correspond to one "bad" and two "good" measurements 
(giving speed-up contributions of and 2 + 2 respectively, totalling 4). For the second measure, 
they correspond to one "good" and two "bad" measurements (giving speed-up contributions of 1 
and 1/2 + 1/2 respectively, totalling 2). 

For a measurement with non-ideal detectors, the dynamics remains isotropic in the Bloch space. 
Moreover, the classification to "good" and "bad" measurements remains the same, as long as the 
detector inefficiency 5 is not greater than 1/2. The inefficiency causes an increase in the purification 
times, that scales as a function of the ratio of inefficiency over impurity, 5/e, in the high-purity 
and high-efficiency limit, (J<1,£<1. Here, the first speed-up quantification, stated in terms of 
attaining a certain average purity, is simply impossible if e < |, and the time required diverges as 
s approaches this limit from above. On the other hand, the second quantification, stated in terms 
of the mean time for individual (stochastically evolving) systems to attain this purity, always gives 
an answer. However we show that the mean time increases exponentially: oc exp [0.5 (5/e] for very 
large 5/e ratio. Still, for moderate 5/e > 1 the mean time is not too much greater than the ideal 
case estimation. 
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